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Abstract 

A model of parton multiple scattering in a dense and expanding medium is described. The simulated results reproduce 
the general features of the data. In particular, in the intermediate trigger momentum region there is a dip-bump 
structure, while at higher trigger momentum the double bumps merge into a central peak. Also, a new measure is 
proposed to quantify the azimuthal distribution with the virtue that it suppresses the statistical fluctuations event- 
by-event, while enhancing the even-structure of the signal. 



1. Introduction 

The study of correlations of jets in heavy-ion col- 
lisions has been a very active area of research for a 
number of reasons. From the theoretical perspec- 
tive back-to-back jets are the most effective means 
to probe a dense medium by jet quenching [1]. In 
experiments the measurement of the azimuthal dis- 
tribution in events triggered by high-p^ particles 
has yielded bountiful data that provide a rich source 
of information on correlations [2], the properties 
of which are largely not satisfactorily explained by 
QCD. It is in that subject where we present two 
pieces of work that are related, but independent. 

First, we present a model of multiple scattering 
of a parton as it traverses a dense and expanding 
medium. An event generator is constructed to re- 
produce the A(j) distribution at intermediate trigger 
momentum so that the recoil partons can give rise to 
deflected jets in addition to the possibility of being 
totally absorbed. It is then extended to higher trig- 
ger momentum to show that the results are in accord 
with the data that exhibit punch-through jets. 

The second part of our presentation addresses the 
issue of how to improve the measure of the event 
structure by eliminating the necessity of making 



background subtraction, a procedure that is thus far 
always used in the analyses of data. We make use of 
factorial moments to suppress the statistical fluctu- 
ations. The proposed method is applied to the anal- 
ysis of the away-side A0 distribution with notable 
success on simulated events. 



2. A multiple scattering model 

In heavy-ion collisions the events with trigger mo- 
mentum in the 2.5 < pr < 4 GeV/c range and the 
away-side azimuthal distribution shown in Fig. 1 [3] 
have stimulated a great deal of attention. There is 
a dip at ip = A(j) — tt = and double bumps at 
tp ~ ±1 rad. This dip-bump structure has been sug- 
gested as a signature of collective response of the 
medium in the form of a Mach Cone structure, initi- 
ated by sonic boom, Cherenkov radiation of gluons, 
or other mechanisms [4,5,6,7]. Events with simulta- 
neous presence of two associated particles approxi- 
mately symmetric about cp = have been observed 
[8] , which support the existence of Mach cone type 
of events. On the other hand, we find it to be plausi- 
ble that there can be a substantial fraction of events 
where the dip-bump structure is due to the effect of 
multiple scattering [9]. 
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Fig. 1. Dip-bump structure in the away side distribution. 
The associated particle angle (p is defined with respect to the 
direction opposite to the trigger. The data points (solid dots) 
are from PHENIX collaboration [3] measured in Au+Au col- 
lision at 0-5% centrality. Momentum ranges for the triggers 
are 2.5-4 GeV/c and for the associated particles are 1-2.5 
GcV/c. 

We assume that each triggered event of a central 
collision is the result of a hard scattering near the 
surface on the trigger side, randomly initiated within 
the region 



£i?o < r < i? , and |A</>| < 



(1) 



where i?o is the radius of the colliding nuclei, £ and 
£i = sin0 c are close to 1. The recoil parton is di- 
rected opposite to the trigger. The hard collision oc- 
curs at i=l followed by scatterings at i=2, 3, 
The medium undergoes Hubblc-likc expansion with 
the transverse radius [10] 

Ri=Roexp[H/(n-l)}, (2) 

and the density is assumed to take on the form 

Pi(run) = ( Ps /Ti)6(Ri-n). (3) 

The central characteristic of our model is that the re- 
coil parton undergoes multiple scattering with ran- 
dom scattering angle chosen from the forward cone 
defined by a Gaussian distribution with a width 

Pi/Ps 

We adopt the BDMPS[11] form for the energy loss 
dE/dx = Kiy/E/E s . (5) 

The stepsize depends on the instantaneous density 
Pi and parton energy Ei . We assume the form 



A., = A s exp(p s - pi)^Ei/E s , 



(6) 



where p s sets the dimensionless scale of the density 
dependence. Eqs(5) and (6) together lead to the re- 
cursion relation for the parton energy 
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Fig. 2. Samlple tracks, (a) Exit tracks, (b) Absorbed tracks. 



E i+ i = Ei[l - nexp(p s - pi)] 2 , 



(7) 



with k = KiA s /2E s . The time interval between scat- 
tering is 



H+l 



Ai/Aj 



If the parton energy is below the cutoff E cut , the 
parton is absorbed by the medium. 

In our simulation, the effective trigger momentum 
is assumed to be at 4.5GeV and the exit parton mo- 
mentum range 0.7 to 2.2 GeV/c for 1 < j>T(assoc) < 
2.5 GeV/c. The shift of 0.3 GeV/c is to account for 
the recombination of exit parton with at least one 
thermal parton which on average carries a momen- 
tum of 0.3 GeV/c. 

Following parameters are fixed at their default 
values: R a = 6f, £ = £i = 0.8. H=0.07 (see ref.[10]), 
E s =5 GcV. _E cu t=0.3 GcV/c. Other parameters are 
adjusted so that the exit track distribution together 
with a thermal background qualitatively reproduce 
the dip-bump feature of the data. The four ad- 
justable parameters are: p s = 0.63, o\ = 0.88 rad, 
A s = 1.9f and k = 0.17. The dashed line in Fig. 
1 is the thermal background due to the energy de- 
posited in the medium [9] . The solid line rises above 
the dashed line due to the contribution from the ip 
distribution of exit partons multiplied by a factor 
of 1.7 to account for the average number of hadrons 
produced from each exit parton through recombi- 
nation of possibly more than one shower parton 
with thermal partons. We regard the agreement 
with data in Fig. 1 as satisfactory. 

Sample tracks are shown in Fig 2. Both plots in 
the figure are superpositions of many events, with 
one track per event. Exit tracks correspond to those 
persistently bending away from the center. These 
tracks are relatively short, and leave the medium at 
ip w 1. The absorbed tracks tend to swing back and 
forth, and are relatively long. It is this underlying 
pattern that leads to the dip-bump structure at the 
intermediate energy. 

At higher trigger momentum the data reveal only 
one central peak as shown in Fig. 3. Since the exit 
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Fig. 3. Comparison between the model fits to the peaks and 
the data. The data points are from STAR collaboration[12] 
in momentum ranges of the triggers: 8-15 GeV/c and of the 
associated particles for (a) 4-6 GeV/c and for (b) > 6 GeV/c. 

partons take more than 10 fm/c to traverse the 
medium, the thermal partons on the far side are 
significantly reduced at such later times. Thus the 
hadronization of the exit partons is dominated by 
fragmentation, when the associated hadron mo- 
mentum is above 4 GeV/c. The predicted hadron 
distributions are shown by the solid curves in Fig. 3. 

Our results indicate that the present multiple 
scattering model can provide a good description of 
the main features of the data, i.e. the dip-bump 
structure at lower trigger momentum, and the 
merging of double bumps into one central peak 
at high trigger momentum. The important conse- 
quence is that only one jet on the away side, if any, 
accompanies a near-side jet, not two. 

3. Use of the background suppressed 
measures to analyze away-side distribution 

So far wc have focused on only the event-averaged 
data. Next wc investigate the implication of the 
event-by-event description of the model. Wc will 
use the factorial moment (FM) approach to analyze 
the data. Consider an event with multiplicity N on 
the away side. We divide the y-region of interest 
into M bins. FM of order q is defined by 

f q = Jj >: J n >'"j - 1) •••(%- g + 1), Uj > q. (9) 
Normalized factorial moment (NFM) is defined by 



F g = f q /fl 



(10) 



The statistical- limit theorem states that if in each 
bin the particle number has a Poisson-like fluctu- 
ation in the large N limit, then F q approaches 1 
for all q, independent of M. The NFM is defined 
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Fig. 4. (a) Illustration of event- by-event fluctuations of F3. 
(b) A comparison of dispersion curves between F2, and F3. 
The i*2 curve is sharper with its peak close to the statistical 
limit of unity. 

for each event, and involves an average over bins. 
In practice N is finite and F q fluctuates from event 
to event. Fig 4(a) displays the fluctuation of F3 for 
some 500 background events. The distribution of F3, 
i.e. dN/dFs vs F3, is shown in Fig. 4 (b) as the solid 
curve. The dashed curve for F2 is sharper and peaks 
closer to the ideal limit of unit. We will see that it is 
the event-averaged NFM that constitutes the basic 
background suppressed measure [13]. 

We now turn to a toy model to illustrate the uses 
of the FM method. Here the signal is defined by a 
cluster of several particles spread over a small ip in- 
terval. We will refer it as a jet. For our model cal- 
culation we use several particles evenly spaced in 
a small ip interval. Three type of events are con- 
sidered. First, the background(bg) events. Each bg- 
event consists of large number of particles randomly 
distributed over the full ip interval of interest. The 
second is background plus one jet (bg+lj) events. 
The jet location is uniformly distributed in a subdo- 
main within the full range. They mimic deflected jet 
plus background events. The third type is (bg+2j) 
events, where the 2 jets are symmetrically produced 
about tp = at small ip region, representing the 
Mach-cone type events. A typical superposed spec- 
tra of all the three types of events are shown in Fig. 5. 

Fig. 6 gives the (F q ) vs M plots. Plot (a) shows 
the average NFM of the background events. Here 
we have (F q ) ~ 1 for the range of M shown. While 
F q varies from event to event, the event-averaged 
quantity (F q ) is approximately the statistical limit. 
Sensitivity of (F q ) to signals can be seen in plots 
(b) and (c). Deviations of (F q ) from unity become 
noticeable for q=3 and 4, especially at high M. 

Factorial moment approach can also be used to 
measure the fluctuations between two (^-regions. 
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Fig. 5. A toy model: Superpositons of ip distributions for 
the production of bg+2j(dash), bg+lj (dash-dot) and bg 
(solid-line) events. 
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Fig. 6. The M-dependence of (F q ) for various orders of q. 
(a) bg, (b) bg+lj, and (c) bg+2j. 

Consider the case where region-I is for f < 0, and 
region-II, for f > 0. Here the difference measures 
the fluctuation. Define 



{D(p,q)) = {\F I q -F I q I \v). 



(11) 



where the pth power amplifies the fluctuations. To 
track the relative normalization one also needs to 
define the corresponding sums 



(S(p, q )) = ((F q T + F q TI r) 



(12) 



The (D) vs (5) plots are shown in Fig. 7. A glance 
at the figure shows that there is a common pattern 
in all plots. More specifically in each plot the bg 
contribution is well localized and suppressed. The 
bg+lj and bg+2j contributions fan out along ap- 
proximately straight lines with bg+lj events having 
larger slopes than those of the bg+2j events. No- 
tice that the coordinate scales increase quickly as 
p and q increase. As mentioned earlier, the bg+2j 
case mimics the Mach cone type of events, while the 
bg+lj case the deflected jet events. One sees that the 
(D) vs (S) plot can be used to distinguish between 
the two types of events without making background 
subtractions. 




Fig. 7. The (D(p,q)) vs (S(p,q)) plots. Rows: p=2,3 and 4, 
columns: q=2, 3 and 4. 

To summarize, we have investigated the use of FM 
method to analyze away-side Acf> distribution. The 
main advantage is that in plots involving FM mea- 
sures signals tend to stand out and the background 
contribution is suppressed. And there is no need to 
make background subtraction. We envision that the 
FM-method has the potential to provide a common 
framework to compare results from different exper- 
iments and various subtraction schemes. 
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